Dual local and non-local cascades in 3D turbulent Beltrami flows 
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We discuss the possibility of dual local and non-local cascades in a 3D turbulent Beltrami flow, 
with inverse energy cascade and direct helicity cascade, by analogy with 2D turbulence. We discuss 
the corresponding energy spectrum in both local and non-local case. Comparison with a high 
Reynolds number turbulent von Karman flow is provided and discussed. 
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Turbulent flows are usually characterized by a wide 
range of scales of motion. Despite their apparent ran- 
domness, these motions organize themselves in a hier- 
archical manner, resulting in power-law spectra for the 
energy as a function of the wavenumbcr k. The index 
of the power law and its range of scale depends on the 
turbulence geometry. For example, in forced isotropic in- 
compressible turbulence, the energy spectrum varies like 
fc -5 / 3 at large scale in 2D and in the so-called inertial 
range in 3D, while it varies like fc -3 at small scales in 2D 
[H]. These features are traditionnally explained through 
local cascades of inviscid global quadratic invariants of 
the Navier-Stokes equations. In 2D, the conserved quan- 
tities are the kinetic energy E = (1/2) < v 2 > and the 
enstrophy fl =< uj 2 > where u = V x v is the vorticity 
and v the velocity. In 3D, the two quadratic invariants 
are the kinetic energy E, and the helicity H =< uj ■ v >, 
that is sign-indefinite, in contrast to energy and enstro- 
phy that are both non-negative. This difference has been 
originally recognized as responsible for the difference be- 
tween 2D and 3D situation @, Q : in 2D, the enstrophy 
invariant effectively blocks the energy cascade from large 
to small scales, resulting in a dual cascade, with energy 
cascading to large scales, and enstrophy to small scales 
[l[. In 3D isotropic incompressible turbulence, a joint 
helicity and energy cascade is rather expected [3j-|5|, in 
which both the energy and the helicity cascade towards 
the small scales. A possible mechanism allowing coex- 
istence of constant energy and helicity flux is the near 
cancellation of the helicity flux between modes of differ- 
ent helical polarity |f|. 

Once the cascade scenario has been clarified, the com- 
putation of the index of the energy spectrum follows 
by simple dimensional argument on the cascade trans- 
fer rate. For example, in 3D, if one assumes that the 
cascade is local and determined only by the energy flux 
e, one gets an energy spectrum E(k) ~ e 2 / 3 fc- 5 / 3 an d 
a helicity spectrum H(k) ~ i5e _1 / 3 fc~ 5 / 3 , where 8 is the 
helicity flux. In 2D, similar assumptions using the en- 
strophy flux r\ instead (when there is no energy flux) 



leads to E{k) ~ r/ 2 / 3 k~ 3 (with logarithmic corrections 
that may result in a steeper spectrum) and an enstrophy 
spectrum Q(fc) ~ rf-l^k^ 1 at the small scales. In con- 
trast, if one assumes a non-local cascade, as suggested in 
at a rate determined by the large scale shear rate S, 
one gets the same spectrum but with a different prefactor 
: E(k) ~ {i]/S)k~ 3 and fi(fc) - (v/S)^ 1 @. Numerical 
simulations rather supports this second scenario Q. 

In the present paper, we consider an intermediate sit- 
uation between the 2D and 3D geometries, in which the 
basic flow is Beltrami, so that vorticity and velocity are 
colincar everywhere and helicity is maximal. In this case, 
the helicity is of constant sign and can possibly block the 
forward energy cascade, resulting in a dual cascade of 
energy and helicity. We then consider both local and 
non-local cases, and compute the resulting energy and 
helicity spectrum. We argue that such a special situa- 
tion is empirically encoutered in large Reynolds number 
von Karman turbulent flows and compare our predictions 
with energy spectra obtained in these flows. 
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FIG. 1: Summary of the two different possibilities for dual 
helicity and energy cascades as a function of the wavemimber 
in a Beltrami flow. Left : local case; right : non-local case. 

Theoretical background and definitions We consider a 
3D Beltrami flow, such that uj = /iv, where \i is constant. 
This flow is an exact solution of 3D Euler equations. It 
also describes the stationary state of a high Reynolds 
number von Karman flow [9[. In such a flow, vorticity 
is maximal, helicity H = J v • oj = J dkH{k) is of con- 
stant sign and we have the equality H(k) = kE{k) at 
the wave number of the Beltrami flow (while in general, 
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we have only the inequality H < kE). This relation is 
the analog of the relation 51 = k 2 E in 2D turbulence. It 
has some important consequences onto the possible cas- 
cade direction. Indeed, if there exists a non-zero flux of 
energy at small scale k > 1, it will be smaller than the 
corresponding helicity flux by a factor 1/k. Therefore, 
at the small scales, the cascade will be mostly governed 
by the helicity transfer rate S and will correspond to a 
direct helicity cascade. At the large scales k < 1, the en- 
ergy flux dominates instead, and we shall have an inverse 
energy cascade. To compute the corresponding spectra, 
we introduce the rate of energy transfer II^; = k ^{k) anc ^ 
the rate of helicity transfer II # = fc 3jffi , where r(k) is 
a characteristic transfer time at wavenumber k. We con- 
sider below two situations. In the local case, r is governed 
only by interaction with wave numbers that are of com- 
parable size. The only characteristic time scale one can 
build in such a case is an eddy turn-over time, given on 
dimensional ground by 

-1/2 

Tloc(k) = I / q 2 E(q)dq I . (1) 



q 2 E(q)dqj 



In the non-local case, transfers between scales are 
dominated by interactions involving two well separated 
scales. This situation is automatically realized for spectra 
steeper than k~ 3 since in that case, the integral defining 
Tioc is dominated by the behaviour near the minimum 
wavenumber, instead of by the behaviour near k. Non- 
local transfers also arise in situations with a well defined 
large scale mean flow, that produces both sweeping and 
shearing of small scales over a characteristic time scale 
S -1 , where S is the large scale shear rate. In such a 
case, the transfer time is given by the time scale of the 
largest scale. This is an intrinsic property of the flow 
which applies to both Energy and Helicity resulting in 

7~non — loc $ • (2) 

We may now proceed to the spectra derivation. Con- 
sider first the small scale energy transfer, for k > 1. In 
such a case, the helicity transfer is dominant. Assuming 
furthermore the existence of an inertial range of scale, in 
which this helicity transfer is equal to a constant <5, we 
get : 



T{k) 



(3) 



In the local case, this leads to the integro-diffcrential 
equation : 

1/2 

= (4) 



k 2 E(k) (J q 2 E{q)dq\ 



with solution 

E loc (k) ~ £ 2 / 3 fc- 7 / 3 , (5) 
(local direct helicity cascade). 



In the non-local case, we get directly |8| : 

E non - loc {k) ~ SS^k- 2 , (6) 
(non — local direct helicity cascade). 

Consider now the large scale transfer. It is dominated by 
energy transfer. Assuming furthermore an inertial range 
of scales with constant energy transfer rate e, we then 
now get : 



kE 

Tl E = -7TT = e. 

T{k) 



(7) 



In the local case, this now leads to the integro-diffcrential 
equation : 



kE{k) ^jf q 2 E(q)dqj 



1/2 



(8) 



(9) 



with solution 

E loc {k) ~ e 2/3 fc~ 5/3 , 

(local inverse energy cascade) . 

In the non-local case, we get directly : 

E non - loc (k) ~ eS^k- 1 , (10) 
(non — local inverse energy cascade). 

These results are summarized in Fig. 1. We see that the 
difference between local and non-local spectrum can be 
tested both via the slope of the spectra, and using the 
prefactor of the spectra, at variance with the 2D situation 
where at small scale, both local and non-local hypotheses 
lead to the same fc -3 spectrum, with different prefactor. 

We now argue that these new energy spectra, resulting 
from a dual energy and helicity cascade, are indeed ob- 
servable in real flows. For this, we consider a turbulent 
von Karman flow, where it has been empirically shown 
that the stationary state is a Beltrami state 0, [l(J . 

Experimental setup We have worked with a specific 
" axisymmetric" configuration: the von Karman flow gen- 
erated by two counter-rotating impellers in a cylindrical 



vessel [111 ]. The cylinder radius and height are respec- 
tively R = 100 mm and H = 180 mm. The impellers con- 
sist of 185 mm diameter disks fitted with sixteen 20 mm 
high curved blades. The impellers are driven by two inde- 
pendent motors at frequencies f± = /2 • The mean forcing 
frequency is / = (/i + / 2 )/2. The flow is divided into two 
toric cells separated by an azimuthal shear layer. In the 
exact counter-rotating regime, in a statistical sens the 
shear layer is located at mid-distance of each impeller. 

In the present paper, we force the flow through disk 
rotation in the (+) direction associated with the convex 
face of the blades going forward [H| . The working fluid is 
either water or glycerol at different dilution rates (respec- 
tively 100% glycerol, 74% glycerol and pure water). The 
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resulting accessible Reynolds numbers (Re — 2irfR 2 v 1 
with v the kinematic viscosity) vary from 10 2 to 10 6 . 

Measurements are done thanks to a stereoscopic Par- 
ticle Image Velocimetry (S-PIV) system. This system 
provides time series of the 3 components of the veloc- 
ity (radial v r (r,z,t), vertical v z (r,z,t) and azimuthal 
v<f,(r, z, t)), and of the azimuthal vorticity component 
Wcp = d z v r — d r v z in the meridian plane on a 58 x 63 
points grid. Time series sampling frequency is in the 
range 1.7 to 15 Hz and includes 2400 to 4200 samples. 
Fluctuations of the velocity and vorticity are then ob- 
tained after substraction of the mean velocity field to 
the instantaneous velocity field. From the correspond- 
ing time series, we compute the 3 dimensions space-time 
power spectrum E(k,2irf), where k is the wavevector 
and / is the frequency. We deduced the spatial veloc- 
ity spectrum E(k = \\k\\) (or equivalently E(l/X), where 
A = 1 / 1 1 k 1 1 is the reduced wavelength) by integrating over 
all directions of k and over all frequencies. Note that our 
procedure defines genuine spatial spectra, independently 
of any Taylor hypothesis. 
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the spectra is observed in the vertical component. This 
is a classical outcome in anisotropic turbulence (see e.g. 



12j ) because of different scaling in the directions paral- 
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FIG. 2: Spatial spectrum E(l/X) in arbitrary units of the 
velocity fluctuations at Reynolds number Re/1000 = 3.2 
(red), 9.5 (yellow), 21 (light green), 41(green), 122 (emer- 
ald), 266 (light blue), 410 (dark blue) and 1022 (purple). 
(Color online): a) Total; b) r-component; c) (^-component; 
d) z-component. The energy spectra have been normalized 
by R 2 f. 

Results The spatial energy spectra for the three com- 
ponents of velocity at different Reynolds numbers are 
shown in Fig. [2] In all cases, the energy is mainly 
concentrated at the largest scales, around A = 0.1 m, 
corresponding to the cylinder radius. In the azimuthal 
and radial components, as well as in the total energy, 
one observes a well defined inertial range at small scales 
A < 0.03 m, and the very beginning of another scal- 
ing range at scales between 0.03 < A < 0.1 m for large 
enough Reynolds numbers. However, no clear scaling for 



lei and perpendicular to the anisotropy. In our case, the 
rotation of the impellers is along the z-direction, defin- 
ing the natural direction of anisotropy. Moreover, we 
have checked that most of the energy is contained in the 
azimuthal velocity component, i.e. in the direction per- 
pendicular to the rotation axis, in agreement with re- 
cent numerical simulations of helical turbulence [12| . In 



the sequel, we therefore focus on the azimuthal veloc- 
ity v'^ and vorticity ui'^ fluctuations, and on the corre- 
sponding energy E<f, =< v£ >= J E(k)dk and helicity 

H=< v y (tl >=jH(k)dk. 

As a first preliminary check, we compare in Fig. [3] the 
azimuthal helicity spectrum, with E ( f ) (l/\)/\. If the flow 
is indeed Beltrami, the two spectra should coincide over 
the range of scale defining the Beltrami flow. We see that 
they indeed look similar at large scale, which is consistent 
with the fact that the mean (large scale) velocity field in 
von Karman is a Beltrami flow [9j. A more quantitative 
test of the Beltrami property would however require the 
three components of vorticity, which is not available in 
our experimental set-up. 

To get a description of the possible cascades, we now 
concentrate on the energy spectra behaviors, through 
the compensated spectrum Q = E$(l/\)/\. At the 
smallest Reynolds number (red curve), Q is maximal 
at A c ps 1/25 m corresponding to the forcing scale. As 
the Reynolds number is increased, we see that for high 
enough Reynolds numbers (Re>~ 40. 10 3 ) Q is develop- 
ping at large scale a constant range starting at A c and 
extending progressively towards the size of the vessel 
A ps 0.1m. This corresponds to a range of fc _1 scaling 
for the energy spectrum, in agreement with the scaling 
expected for the inverse cascade in the non-local case. 
Below A c , a crossover for Q is observed at any Reynolds 
number. This corresponds to the transition to the direct 
cascade regime where a power law can be measured. We 
observe that spectra at small scales are gradually filled 
with increasing Reynolds number from Aw 1/50 m up to 
reach the resolution limit at A ps 1/130 m. 

The scaling slope in the direct cascade regime is com- 
puted using linear fit with adaptative boundaries. The 
resulting slopes are shown in Fig. |4] for a large range of 
Reynolds number. 

The first striking result is that for any Reynolds num- 
ber, the direct cascade slope is steeper than the clas- 
sical Kolmogorov -5/3 expected in the homogeneous 
and isotropic 3D turbulence. At small Reynolds num- 
ber (Re < 50. 10 3 ), the slope is nearly —2 close to 
the non-local prediction. The slope then steepens fur- 
ther at increasing Reynolds number reaching —2.2 at 
Re > 100.10 3 , close to the local prediction (-7/3). At 
the present time, our time statistics (3600 frames) does 
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FIG. 3: Top: azimuthal contribution of helicity H — j v^uj^. 
Bottom: spatial spectrum Q = E t / > (l/\)/\ of the azimuthal 
velocity fluctuations. The spectra have been arbitrarily 
shifted for clarity. Same color code as in Fig. [5] Lines are lin- 
ear fittings in the direct cascade range 25 < A -1 < 50 (resp. 
30 < A -1 < 100) when Re < 41000 (resp. Re > 65000). 

not allow us to restrict the error bars on the fit below 10 
percent. This is enough to exclude the Kolmogorov case 
(exponent —5/3) from our results, but not to assess with- 
out ambiguity the observed shift from a " non-local" to a 
"local" spectrum at Re ~ 5.10 4 . It is however interest- 
ing to note that the transition Reynolds number is pre- 
cisely associated with a spontaneous symmetry breaking 
bifurcation of the mean flow [ll[, with possible associ- 
ated changes in the structure of the vortices of the shear 
layer. Work is in progress to check whether this is likely 
to change the mean shear felt by the smallest scale, and, 
therefore, the locality properties of the cascade. 

Discussion We have studied spatial energy spectra in 
a von Karman experiment, in a wide range of Reynolds 
numbers. In all cases, we have observed a dual en- 
ergy cascade, one at the scales larger than the forcing 
scale, with increasing energy concentration at increas- 
ing Reynolds number, and one at the scales smaller than 
the forcing scale, with increasingly steeper slope as the 
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FIG. 4: Slopes of E(k) for various Reynolds numbers (color 
online). Green crosses, black stars and blue circles corre- 
sponds to resp. 100%, 74% glycerol mixture and pure water. 
Error bars are computed using the same time series as in 
Fig. [21 Horizontal black lines mark the local (-7/3) and non- 
local (-2) predictions. The dashed line marks the Kolmogorov 
(-5/3) index. 

Reynolds number increases. None of these cascades scale 
according to the fc -5 / 3 Kolmogorov prediction, valid for 
isotropic homogeneous turbulence. On the other hand, 
their slope is in agreement with the slopes predicted using 
the helicity conservation in a Beltrami like flow. Given 
that the von Karman experiment is neither isotropic, nor 
homogeneous, the discrepancy with Kolmogorov predic- 
tion is not surprising. In addition, its Beltrami character 
at large scale [§] makes it suitable to test the theoreti- 
cal prediction. It is nevertheless interesting to compare 
our experimental findings with other non-isotropic sit- 
uations, where non-Kolmogorov spectra have been pre- 
dicted and observed. For example, in the rotating he- 
lical case, Mininni and Pouquet 12[ observe a similar 
coexistence of a large scale inverse energy cascade, and 
a direct cascade towards small scales. The slope of the 
energy and helicity spectra in the direct cascade is close 
to —2, like in the present experiment. However, using 
argument adapted from Iroshnikov-Kraichnan [13| . they 
predict that for flow with maximal helicity, the scaling 
should be k~ 5 ^ 2 for the energy, and fc~ 3 / 2 for the he- 
licity. This is not observed in our experiment, although 
we are in a situation with large (but perhaps nor com- 
pletely maximum) helicity. It would be interesting in 
that respect to perform new von Karman experiments 
with added rotation, to see whether a steepening of the 
energy spectrum occurs. 
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